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Abstract 

We derive general expressions for the sum rules of the eigenvalues of drums of 
arbitrary shape and arbitrary density, obeying different boundary conditions. 
The formulas that we present are a generalization of the analogous formulas for 
one dimensional inhomogeneous systems that we have obtained in a previous 
paper. We also discuss the extension of these formulas to higher dimensions. 
We show that in the special case of a density depending only on one variable 
the sum rules of any integer order can be expressed in terms of a single series. 
As an application of our result we derive exact sum rules for the a homoge- 
neous circular annulus with different boundary conditions, for a homogeneous 
circular sector and for a radially inhomogeneous circular annulus with Dirichlet 
boundary conditions. 
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1. Introduction 

This paper considers the problem of obtaining explicit expressions for the 
sum rules 

^^P^ = Y.^ ' P = 2,3,... (1) 

n 

where En are the eigenvalues of the Helmholtz equation 

{-A)ipn{x, y) = £^„S(x, y)'ipn{x, y) . (2) 

over a two dimensional domain (x,?/) € fi. We assume that 17 is a domain 
where a orthonormal basis is known (square, circle, etc) since the more general 
problem on an arbitrary domain may be reduced to the form of eq.Q using a 
confornial map. Here n is the set of quantum numbers which fully specifies a 
solution and Y,{x,y) may either be a physical density, a density obtained from 
a conformal map or a composition of the two. 
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Unfortunately, the evaluation of Z{p) using eq.dT]) requires the knowledge 
of the eigenvalues of the problem, which is granted only in special cases. In a 
recent paper, ref.[H, we have discussed an alternative approach which uses the 
representation of Z{p) directly as the trace of the hermitian operator 

Qp = (Vs(-A)-iVs)'' (3) 

where p > 1 is a real exponent. We have proved there that the eigenvalues of 
Q are the reciprocals of the eigenvalues of eq.([2]) and therefore the trace of 
provides the spectral zeta function associated with the problem of eq. ([2| . For 
integer values of p it is possible to obtain explicit expressions for the trace and 
therefore for Z(p): we have discussed specific examples in ref. In the general 
case of non integer values of p, we have shown in ref.[l| that it is possible to 
evaluate the trace using perturbation theory; although the expressions obtained 
following this approach are defined for p > 1 (in d dimensions, p > (i/2) it 
is possible to perform their analytic continuation to s < 1: in this way we 
have obtained the Casimir energy of specific systems in one, two and three 
dimensions. 

There is a large number of works where spectral sum rules for given sytems 
have been studied: these include the sum rules for quantum mechanical an- 
harmonic oscillators 0, Si, for potentials V{t) = gr^ (p > 0, g > 0) jik for 
Aharonov-Bohm quantum billiards Q , for the zeroes of Bessel functions 4 J3 j 
for the Selberg's zeta function for compact Riemann surfaces Q (see also J), 
for quantum mechanical one dimensional potentials [lol|. for two dimensional do- 



13|,ll4|. 



mains close to the unit disk 11 1, for the cardioid and related domains jl2|. for a 
T^T-symmetric hamiltonian 



In this paper we provide a method to calculate the sum rules for the eigen- 
values of the negative laplacian in two and higher dimensions, for systems of ar- 
bitrary density. The approach that we describe is based on the results of ref. [if , 
which are now expressed in terms of the appropriate Green's functions and used 
to derive a general formula for the sum rule of arbitrary order p > d/2 [d are 
the dimensions of the problem) . This formula is a generalization of the analo- 
gous formula for one-dimensional inhomogeneous systems obtained in ref. (isj . 
where it has been applied to a number of examples (in the case of a particular 
inhomogeneous string we have been able to obtain sum rules up to order nine). 
In two dimensions the reference domain chosen for our calculation is a rectangle 
(for d > 2 we use its generalization), although it can be any domain where 
a Green's function is available. In particular in two dimensions, our formula 
should be compared with the equivalent formula obtained by Itzykson, Moussa 
and Luck [l^ , who expressed the sum rule of a given order n for the Dirichlet 
eigenvalues of the negative laplacian on an arbitrary domain T in terms of the 
Green's function on the upper half plane; moreover Dittmar [l3] has obtained a 
formula for the sum rule of order two for the eigenvalues of the negative lapla- 
cian on a two dimensional domain both for Dirichlet and Neumann boundary 
conditions in terms of the Green's functions on a circle (see Theorems 2.1 and 
3.3 of ref. 17]). The expressions that we derive in this paper apply to sum rules 
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of arbitrary order, to different boundary conditions an both to homogenous and 
inhomogeneous systems. 

The paper is organized as follows: in Section [5] we derive the expressions for 
the Green's functions on a rectangle of sides a and 6, with different boundary 
conditions; in Section [3] we obtain a general formula for the sum rules of arbi- 
trary order in two dimensions; in Section 2] we discuss the generalization of this 
formula to higher dimensions; in Section[5]we consider some applications of our 
formula; finally, in Section [5] we draw our conclusions. 



2. Green's functions on a rectangle 

In this section we derive explicit expressions for the Green's functions of the 
negative laplacian on a rectangle of sides a and b and obeying different boundary 
conditions at the borders. We consider the cases of Dirichlet, Neumann and 
periodic be on the borders of the rectangle, of mixed be, Dirichlet-Neumann, 
Dirichlet-periodic, Neumann-periodic and mixed Neumann-Dirichlet boundary 
condition in one direction and periodic in the other direction. 

2.1. Dirichlet boundary conditions 

We consider the rectangular region {—a/2, a/2) x (— &/2, 6/2) and assume 
Dirichlet boundary conditions at the border of this region; the Green's function 
for the negative laplacian is 

n^ = lny=l ^ri^ + Vriy 

Here V'n^ (x) and </)„^ (y) are the Dirichlet eigenfunction of the ID negative 
laplacian on x £ {—a/2, a/2) and y e (—6/2,6/2) respectively: 

i'n^ix)= sm(J^{x + a/2)^ ' ^ 1.2,... (5) 

V'„Jy)^y|sin(^(y + 6/2)) , Uy ^ 1,2,... (6) 

ei^'' and ryi^^* are the Dirichlet eigenvalues of the ID negative laplacian in 
the two directions: 

AD) ^ !^ <D) ^ ^ 

"x ^2 ' '"y 6^ ' ^ 

Eq.(j4]) may be cast in the form: 

oo 

G^°\x,y-x',y') = ^ gif ^2/, yOV-i^^^ WV-i^^H^') , (8) 
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where (see Eq.(3.169) and problem 2.15 of Rcf. 

i/n^ \y^y ) — 2^ (£1) 



n„ = l 



sinh [\[lS}{y^ + sinh [\[W}{hl2 - y>) 



ei^^sinh\/ei^-'6 



(9) 



We have introduced the notation ?/< — min(y,j/') and j/> = inax(y,?/'). 
We also define 

G(^)(a;,y;x',y') ^ ^ (a;, y; x', - 2/') 

+ GL''^(a;,y;x',y')^(y'-y) , (10) 

where the explicit form of G[^\x,y]x\y') is easily obtained from Eq. (|S]). 
Notice also that 

G^^\x,y-x',y') = Gi^\x,-y;x\-y') 
Gi^\x,y';x',y) = Gi^\x,y; x' ,y') . 

It is also possible to express Eq.® in an alternative form, observing that 

^ =2Y^ p-(2i+l)7™6/o (II) 

sinh HH^ ^ ' ^ ' 

a j=0 

Upon substitution of this expansion inside Eq.® we have 

oc 

^]{x-,x+,y-,y+) 



where 

, ^ f , f T^C^tjj + b~ y+) \ ('. 
ilj[x^,x+,y^,y+) = I cosh I I — cos 



' ^{2h]+b - y+) 
V V a 



l^cosh ^ 

cosh (^i^^^l±ll-MU cos r^-^ 



/nx 
\ a 

/TTX 



4 



cosh 



7T{\y-\+2bj) 



a 



( cosh (^i^^^A±i^±)]+ COS (^) 

cos 

V a 



and x±=xi± X2, y± = yi± y2- 
Notice that 

eo(0,a;+,0,y+) =0 

and therefore G^^^ (x, y; x' , y') diverges logarithmically when X2 — > Xi and j/2 — > 
2/1- 

Although the series in Ea. (jl2p converges more rapidly than the series in 
Eq.®, in general the latter is more appropriate for the evaluation of the sum 
rules. 

2.2. Neumann boundary conditions 

We now come to the calculation of the Green's function for Neumann bound- 
ary conditions 0; in this case we have 

u [x,y,x ,y ) — jj^^ [.^0) 



rij: ,Uj- .riy .Utj 



where u n « is the sum over all possible values of the quantum number, 
with the exclusion of the divergent term corresponding to n^; = = and 
Ux = Uy ^ 1 H 

Here ei^lu^ and 7}n^]uy are the Neumann eigenvalues of the negative laplacian 
in the two orthogonal directions; 

2 _2 

Ux = 2 , nx = l,2 



4!k - \ (.37-.- ".-0,1,2,... (^g^ 



(N) ^ J — p— ^ , My = 1 , nj^=0,l,2, 

, '^y — 2 , Hy 1,2,. 



''^ny,Uy S {2ny~lfTr^ 



^We have not found references reporting the equivalent expressions of eq. ||9ll for boundary 
conditions other than Dirichlet. The results that we report for these cases have been derived 
in the present paper following the same approach used for eq. lO. 

^ Erix, tlx, "h, '^""'•"^ '"«•"!' = Enx=l J^iix = l '^"•^'""I'O'l ^"h = 1 ^i1h = 1 '^'''^■"!/>"h 

nx— 1 ^— JUx — 1 ^-^riy—l Z_/Uy — 1 ^X jt^X it^y ,Tiy • 
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and '0n^,Mx i^) and (jjO^My (j/) are its eigenfunctions 



2 cos^^i^ 

a a 



2 ■ (2n^—l)Trx 
sm 



2 2nT,7r'V 



1 „• (2Tt„-l)7ry 
b ''"^ fc 



Tlj; = , Ux = l 

> , = 1 (17) 

TIj: > , Ux —2 

Uy = , Uy = 1 

Hy >0 , Uy ^ 1 (18) 

rij^ > , Uy = 2 



In this case Eq. ([T5t may be expressed as 



GW(x,2/;x',y') 



oo 2 

E E & iy^ (^XL (^') (19) 



riT — l Ut — I 



where gif^'ju^ {y, y') is the analogous of Eq.(|9]) for Neumann boundary conditions: 

l-6„ 



a^ZlSv^y') 



n^.o , ^ ^ (?!>i^,l„(y)'^i^,l„(y') 



, W ^ ^ (AT) (AT) 

Utn^,u^ Uy = lny = l <^n^,u^ ^ ilny,Uy 



For rij; = and itj; = 1, g'nP,u.^{yTy') reduces to the one-dimensional Green's 
function with Neumann boundary conditions, which is reported in Ref. il5]: 



9iV{y.y') 



9 + Qb{y<-y>)+Q{yl+yl) 



126 



(20) 



For Ux > Q and u^ — 1,2 we have 



Uy,v') = 



cosh ( V£L^ ) cosh ( Jv'JlL^ 



r^nju, sinh ( J T]i^}uJ 



(21) 



2.3. Dirichlet- Neumann boundary conditions 

We now come to the calculation of the Green's function for mixed Dirichlet- 
Neumann boundary conditions; in this case we have 



G^^''\x,y;x',y') 



I'^y 



^(D) (N) 



6 



where 

oo c« oo 2 



(^n^,ny,Uy — ^ C„^,o,l + ^ ^ ^ C„^,„j,,„j, . (23) 



1 riy — 1 ^iy— 1 



We may express eq. (|22l) as: 

-.{DN} 

~b 



{DN}i M oo 2 



riy — 1 Uy— 1 

where 

For riy = and = 1 gr!^^l{x, x') is the ID Green's function for Dirichlet 
be reported in Ref.[15|: 

g\^f\:cy) = («-2^>)(« + 2^<) (26) 
while for ny > and Uy — 1,2 

sinh ( J 'q'^njuy (x> - f ) ) sinh ( J Vn^yi^ + x<) 

for Uy > and Uy = 1, 2. 

Alternatively we may write the Green's function as 



G(^^)(x,2/;x',y') = (28) 

where 



{y)4>ny.uy {y') 

tiy = lny = l *^nx ~r llny,Uy 

cosh (J^\b/2 - yy)) cosh (J^{b/2 + y<) 



el^) sinh(feVe(^)) 



(29) 
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2.4- Periodic boundary conditions 

We now come to the calculation of the Green's function for periodic boundary 
conditions; in this case we have 

iyWnJuAx')<t>ny:Uy{y') 



'^y 



IP) , 

^nx,Ux ' 'lny,Uy 



where we have excluded from the sum the divergent term corresponding to 

rix — riy — Q and = Uy = \ as for the case of Neumann be. 

(p) (p) 

Here e-nj.u^ and r/uy 'uy are the Neumann eigenvalues of the negative laplacian 
in the two orthogonal directions; 



and ipnj,u:^{x) and (f>ny,uy{y) are its eigenf unctions: 



AP) 



2^2 



V 



62 



(31) 



2 sin 



2cos^ 

2sinHnp^ 

b 



Ur, > , Ux 

nx>0, Ux 



Uy — , Uy — I 
Uy > , Uy = 1 
Uy > , Uy =2 



(32) 



(33) 



In this case Eq. (PU)) may be expressed as 



where 



rij. — 1 U3. — I 



1 - 5n^.fl 



^ — (t>i^y]uyiy)(f'lry]uyiy') 



E E 



AP) 



AP) 



For Ux — and u^ — I, sl^L^l?/,?/') reduces to the one-dimensional Green's 
function with periodic boundary conditions reported in ref. [l5| : 



b^-6b\y-y'\+6{y-y'f 



12b 



For Ux > and Ux — 1,2 we have 



cosh 



v/^(|2/-y'|-6/2) 



sinh 



(35) 



(36) 
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2.5. Dirichlet- Periodic boundary conditions 

We now come to the calculation of the Green's function for mixed Dirichlet- 
periodic boundary conditions; in this case we have 

'''''' '--^-''luMj^-"'^'''^ 

(D) _^(P) 



riy ,Uy 



where V„ „ „ is defined in eg.d^. 

We may express this Green's functions as 

>{x,y;x,y) = - 



b 

cxj 2 

E E (y ) (38) 

riy—l Uy—1 



where 

JDP 



9^Z^M - E^^#4^ (39) 



For Tij, = and Mj, — 1, 5n^,u), (x, x') reduces to the ID Green's function for 
Dirichlet be: 

g^^P = («-2x>)(a + 2.<) ^^^^ 

For Tij, > and Uj^ = 1 , 2 

sinh f y5S^(x> - f )) sinh (^^y{^ + .t<)) 

<^^(x,x') = ^ ^^''^ 

V^yi^k sinh Uiii^yluy a\ 

Alternatively we may write the Green's function as 

oo 

G(^^)(x,2/;x',y') = E CHy'2^')V^"^^(^)CH^') (42) 

where 



2 oo 

^rH.,y') - -4i5T+E E 



^ °° (t}^fy]uyiy)(l>';fy]uy{y') 

{D) ' ^ ID) . JP) 



coshfv^(|y-2/'|-V2) 



(43) 



2\/ei^^ sinh ( 
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2.6. Neumann-Periodic boundary conditions 

We now come to the calculation of the Green's function for mixed Neumann- 
periodic boundary conditions: 

where X^n « n « been defined earlier for the case of Neumann boundary 
conditions. 

We may express the Green's functions as 

(NP)f ,^ oo 2 



where 

yny,Uy\-^^-^ ) — (P) ^ (N) (P) • 

1 11ny,Uy nx = ltix = l ^ri^,u^ ~r rjny.Uy 

For nj, = and My = 1, .gi^.u! (^i reduces to the ID Green's function for 
Neumann be reported in Ref. [15j 

5^7) (a:, a:') = - ^^\^ - -'^+ ^i-' + ^'") _ (47) 

For fly > and Uy — 1,2 



cosh ( v/.i:k ^ cosh J.i^k ^ 



C{;:(2^,2;') = ^ ^ ^.(48) 

Alternatively we may write the Green's function as 
G^^^Hx,y,x\y') = ^oi (^'^) + ^ E ^^^2 (^/^ ^OCL (^)<U(-^) 



where 



-, c 2 00 ,(P) I -sAP) I l\ 

n'-^PUii - i_lf!!i£ , \^ (Pny,Uy(y)(pny,Uy(y ) , , 

9n^,uSy^y > - (N) +2^ 2^ (jv) , (P) • 

erix.ttx Uj,=lnj,=0 enx,«x + ^"y,"!/ 

For riy = and My = 1, q ^.m ! (x, reduces to the ID Green's function for 
periodic be reported in Ref. [l5| 

9Z'\y^v') - (51) 
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For rij, > and Uy — 1,2 



cosh fy^(|y-2/'|- 6/2)) 

= ^ ^ 



2Vei-) sinh ^ 



2.7. Neumann-Dirichlet- Periodic boundary conditions 

We now come to the calculation ol the Green's function for mixed Neumann- 
Dirichlet-periodic boundary conditions: in this case we assume Neumann be at 
X = —Objl, Dirichlet be at x = +a/2 and periodic boundary conditions at 

y = ±6/2. 
We have 

r<{NDP)i I l\ _ V'wx ^ {^)'t>'ny,Uy(y)i'n^ \x')(ikx^.Uy{v'), ^^^ 



i^y I'^y 



^71^ ~r ^7ly ,Uy 



where V]^ ^ „ is defined in eq. 
We have 



and 



C^H.) = V!-- ( "^'"^"If"^'"^ ) (53) 



4a2 

We may cast the Green's function of Eq.(|52|) as 

(NDP), /X oo 2 

G(^^^)(x,2/;.',y') = V ^ ^ <.gr(x,xO<!.„(y))C4(^) 

rly — 1 Uj^— 1 

where 



For rij, = and Mj, = 1 we have the one-dimensional Green's function for 
mixed Neumann-Dirichlet boundary conditions: 

5^,;i^^'(x,x') = {-x^+a/2) (57) 
For Uy > and = 1 , 2 we have 

sech ( aJelZL ] sinh 



\/erC!«„^ sinh {\\f(^^y{a - 2a;>)^ cosh ^^\/el^!«Ja -f 2a;<) 



cr(-'-') = — ^ ' — ^ — ^ ^ 

^ny,Uy 
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Alternatively we may cast the Green's function of EQ. f[5^ as 

oo 

G(^^^)(x,y;x',y') = E C'^Hy, 2/')^^^))^^-') (59) 



where 



Un^ (yTy ) — ^ ^ND-j 2^ 2^ {ND) (P) 



^ ^"X ' ny=lUy=l ' Vriyl 

cosh(^V^(|y-y'|-V2)) 



ND) 



(60) 



2.8. Dirichlet- Neumann- Periodic boundary conditions 

In this case we assume Dirichlet be at a; = —a/2, Neumann be at a; = +a/2 

and periodic boundary conditions at y = ±b/2. This case is trivial since the 
basis with 

^i^^Hx) = C^'H-^) (61) 

Therefore 

ffiX^(x,x') = gi^Jpi-^,-:c') (62) 

3. Exact sum rules 

In Ref. [13] we have derived a set of rules which allow one to obtain explicit 
expressions for the sum rules for the eigenvalues of an inhomogcneous string: 



E 4 ' (63) 



p p 

with n = 1, 2, 

These rules may be easily generalized to higher dimensions: in particular in 
two dimensions the sum rule of order n may be expressed as 

4^' = j d^Ri J d^R2 ...J d^Rn G(R„; Ri) 

n — 1 n 

• n [^(^-^'^+i)]np(^*)] (64) 

where Ri = {xi,yi) a.nd\xi\ < a/2, \yi\ < b/2. Here G'(Ri; Rj) is the 2D Green's 
function obeying the appropriate be and S(Ri) is a density 



•^Notice that G(R, R') is the Green's function on any domain where a basis is known, such 
as a rectangle or a circle. 
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As done in Ref. 15| it is convenient to cast this expression in an alternative 



form, using the "y-ordered" Green's function: 

/a/2 /■a/2 /■a/2 /■b/2 /■yi rVn-l 

dxi j dx2... dxn / dyi dy2 ■ . ■ dy^ 

-a/2 J-a/2 J -a/2 J -b/2 J-b/2 J-b/2 

n 

■ g(Ri,R2,...,R„)nP(R0] > (65) 

i=l 

where 

e(Ri,R2,...,R„)= |g(R„;Ri)[] [G(R,; R,+i)]| (66) 

I i=l ) -p 

and 

{/(Ri,...,R„)}p= J2 /(Rpi,---,RpJ • (67) 

permutations 

For example, the expressions for Q up to order 4 are: 

g{Ri) = G+{xi,yi;xi,yi) 

e(Ri,R2) = 2 [G+(xi,yi;x2,y2)]' 

fJ(Ri,R2,R3) = 6 G+{xi,yi;x2,y2)G+{x2,y2;x3,y3)G+{xi,yi;x3,y3) 

C/ (Ri,R2,R3,R4) = 8 [G+{xi,yi;x2,y2)G+{xi,yi;x4,y4)G+{x2,y2;x3,y3)G+{x3,y3;x4,y4) 

+ G+{xi,yi; X3, y3)G+{xi,yi; x^, y4)G+{x2,y2; xs, y3)G+{x2,y2; xa, 2/4) 

+ G+{xi,yi; X2,y2)G+{xi,yi;x3,y3)G+{x2,y2; X4, y4)G+{x3, ys; X4, 7/4)] 

We can therefore calculate Z{n) with the diagrammatic rules: 

• Draw n points Ri , . . . , Rn on a line; 

• Connect each point to any two other points in all possible inequivalent 
ways excluding the disconnected diagrams and the diagrams corresponding 
to a cyclic permutation of the points; 



Associate a density S(Ri) at each point R; (i = 1, . . . , n) 



• Associate a factor G+(Ri; Rj) to each line connecting to Rj (i < j); 

• Multiply the result by a factor 2n, corresponding to the n cyclic permu- 
tations of each inequivalent configuration and to the 2 possible directions 
in which each diagram can be traveled; 

• Integrate the expression obtained from the steps above over the internal 
points: 

/■a/2 /■a/2 /•a/2 /■b/2 /■yi ry-n-l 

/ dxi f dx2... I dxn 1 dyi dy2 . ■ . dyn 



a/2 J-a/2 J-a/2 J-b/2 J-b/2 J-b/2 
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It is easy to convince oneself that for the sum rule of order n there are 
n!/2n = (n — l)!/2 independent diagrams (for n > 2). 

These rules imply that the sum rule of order n for a general density contains 
n series, one for each factor G'+(Rj;Rj) appearing in the expression; however, 
in the case of a density which depends only on one variable the orthogonality of 
the eigenfunctions along the homogeneous direction allow to reduce the multiple 
series to a single series, leaving only the integrals along the direction where the 
density varies. 

Assuming for simplicity that S = S(y), and expressing for the Green's func- 
tion as 

oo 

G+{x,y;x',y') = ^5^+)(2/, 2/')V'p(a;)V'„(a;') (68) 



/b/2 nyi rVn-l " 

dyi / dy2... dyn Q{yi,y2,...,yn)l[ [S(2/i()$9,) 
-6/2 J-b/2 J-b/2 ~.i 



we obtain: 

fb/2 pyi rVn-l 

dyi / dy2... 

p ^-b/2 J-b/2 J-b/2 

where 

Q{yi,y2,---,yn) = <gp{yn;yi)Yl[gp{yi;yi+i)]> ■ (70) 



4. Higher dimensions 

We briefly discuss how these results generalize to the case of > 2 dimen- 
sions; we consider a d-dimensional region with \xi\ < ai/2 and i = 1,. . . ,d. For 
simplicity we restrict our analysis to Dirichlet boundary conditions, since the 
cases corresponding to the other boundary conditions can be obtained in an 
analogous way. 

The Green's function reads 



G^'^\xi,...,Xd;x[,...,x'^) - (l-D) , , (dD) 

ni = l ntj = l £"1 

• 4\^)(a.r) . . . Vr)(x.)4\^)(xi) . . . Vr)(23l) 

where: 



and 



^C^^)(x,)EE,/^ sinf^(x. + a,/2)') , = 1,2,. 



e<i^^^^ , n, = 1,2,... . (72) 

CI A 
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We may now apply Eq.® on any of the d series contained in Ea. ([7T|) : for 
instance, if we use it on the d*'' series we have 



oo oo 



"1 = 1 «d-l = l 



where 

„(D) / ^ V \ ^ \^ ^"<i yXdjVn^ [Xj^} 

yni,...,na-i\-^<i,Xci) — {dD) ^ 



sinh (yFni,...,„^_i(xd< + ad/2)) sinh (^r„^^...^„^_^ (0^/2 - '^d-^^^ 

•\/rni,...,Tid_i sinh n\,...,nd-x'^d 



and 



1=1 

The expressions for the sum rules in this case are the analogous of the ones 
discussed in the previous section, with the appropriate Green's functions and 
density, which is now a function of d variables. 

It is interesting to see what happens in the case of a density which depends 
only on one direction, for instance S = iMx^. In this case, as for the two dimen- 
sional case, the eigenfunctions corresponding to the homogeneous directions can 
be eliminated from the expressions for the sum rules, using their orthogonality. 
Therefore the sum rules of any integer order for a d-dimensional system of this 
kind can be expressed in term of d — 1 infinite series. 



5. Applications 

We now discuss few applications of the results obtained in the previous 
sections. 

5.1. Circular annulus 

In Refs.P, we have discussed the application of perturbation theory to 
the study of the spectrum of a circular annulus. In particular, in ref. [ij we have 
evaluated the sum rule of order two for the eigenvalues of the circular annulus 
of unit external radius and internal radius r, using the matrix elements of the 
conformal density and we have calculated the Casimir energy of the annulus in 
the limit r ^ 1^. 

We will now apply the results of the previous sections to obtain explicit 
expressions for the sum rules for a circular annulus: the function 

/(z) ^pZ+liog-™- (76) 
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Figure 1: Circular annulus obtained acting with the map I I76I I the rectangle 
log r-min, —5 log i"min] X [^T, 7r] . Periodic boundary conditions are imposed on the dashed 
line. 



maps the rectangle \ogrmin, — \ logr-,„i„] x [— tt, tt] onto a circular annulus of 
external radius R = 1 and internal radius rmin ■ Instead of solving the Helmholtz 
equation on the annulus, using the conformal map we may solve an equivalent 
Helmholtz equation for an inhomogeneous medium of density 

E(a;, y) = r™,„ 6^"= (77) 

over the rectangle [l|,[l^[23l- Since this density depends only on one coordinate, 
we will be able to obtain general expressions for the sum rules of the circular 
annulus in terms of a single series, as we have pointed out in the previous section. 

We will discuss the two cases of an annulus with either Dirichlet or Neumann 
be at the border: the first case corresponds to using Dirichlet be along the x 
direction and periodic be along the y direction; the second case corresponds to 
using Neumann be along the x direction and periodic be along the y direction. 
Notice that using Dirichlet be alongthe y direction would correspond to studying 
an annulus with a transverse cut [l|. 

Applying the general formulas derived in this paper it is possible to obtain 
explicit expressions for the sum rules of the circular annulus; for example the 
sum rule of order two is obtained as 



- log rrr 



Z2 {r,nin) = 2 / dx dx' 5q ;i {x,x') J^{x)T.{x') 



2 



16 



^ log Tmin r- -i 2 



Using the explicit expression for g^Y"^ {x, x') and performing the integrations 
one obtains 

Z^^^Vr ■ ) = ^m.n log' (^min) _ 1 _ ^^ .2, . ) , 1 _ if / . 



^i^n log'(r^in) _ log^(r^m) .^g^g ^ ^3^6^ ^ g2^; 

2(^L„-1) 144(rL„ + l) 
+ 73r^„ - Stt^ (r^,„ + l)^ (r^,„ + l) + 26) 

n (n2 + 5) (r^,„ - l) r^",„ (^2^^^ - l) 



6 _L «0».4 

min 



n— 3 



min 



; („2 _ 4) („2 _ 1)2 (^2^„^ _ 1)2 3 _ 1)2 (^2„^^ _ i)2 



(78) 



This expression may be cast in terms of a more rapidly convergent series 
observing that < Tmin < 1 and expanding the factor 1/ (r^^„ — l) in powers 
of r. After performing the summation over n one is left with the new series: 

7iDP),^ . ^ _ 
■^2 V' mm } — 



(''mm ^ l) ('nwn + l) ^12 (''mm) . (''mm ^ l) (^min + l) ^12 (^fnm) 



16r^ ■ 8r^ ■ 

\2 



mm 



log(r„,„,) 5 (r^„ - l) (r^i„ - l)^ 



4(l-?-mm) 641og(r™„) 161og2(r„i„) 
{rLn - 1)' (rLn + 1) log (1 - 



8r^ . 

min 



+ 



16r-lm 2 (r2„„ - 1) 



1 



"I" gyg ['''l^min 202r^j„ + 274r^j„ 297r^j„ + 127r 



2/4 

'r, 



mzn 



36 «« 2 ^ 4 (5r^„ - 54rl,„ - 27) 



+ 66r„,,„ + 



^min (''mm + 1) 



°° r 1 

' / j \J ' / ig V mm / y^' mm y mm y mm J 

_ 1(,.4. - 1) - lV-'^^+'^Li2 fr^^+^'l 

+ 1(^2 _ 1) V^4j+4 ^\ -20+1) / 2j+2\ 
' g V mjra -^z mm ' J ' mm ^' mm J 

+ ^ fr'' - 1) (r^^- - fr^^- + r^^ - + l^l r"^^ log f 1 - r^^ 



16 

— fr^ . - 1) ir''+^ - iV ir^'+* + r"^+' + l") r-^^^^'^ log fl - r'^+^'l 

Ig V mm ^ J y mm J y mm ' mm ' J mm 6 \^ mm J 
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1 



V mm 

-20+2) 

576 



(-3(rL„-l)'( 



-2(j+l) 



log (l 



mm I 



4r 



1 r 



ij+i 



+ 36(r^,„-l)') 



Clearly the rate of convergence of this series increases for rmin 

^ 0+: in this 

limit one obtains the behavior 



48 



5_ 
32 



1 



16 log (rmin) 

7 



64 log (r,„i„) 



+ 



81og^(r„i„) 48 



Notice that 



lim z'^^^Ur ■ ) 

■ ^2 \' mm) — 



TT 

48 



5 

32 



0.04936675836 



(80) 



(81) 



is the Dirichlet sum rule of order two for the unit circle. 

Notice also that the derivative of z'^^\rmin) diverges at r = 0: 



lim 

rmin->-0+ 



1 



S^min log (ryTiiyi) 64r777,i7T, log (v^mm) 



= -oo (82) 



In a similar way we may obtain the sum rule of higher orders: although we 
have proved that these sum rules will involve a single series, since the density 
depends only on one direction, their expressions become lengthier and we do 
not see any advantage in reporting them here. 

Instead, we investigate explicitly the limit of an infinitesimal hole, as done 
in the previous case: 



3 \' mm) 



C(3) 3^ 
32 768 



TT 

128 



+ 



15 



23 



64 log^ {rmin ) 512 log^ {rmin ) 
3 15 

+ 



11521og(rTOi„) 
19 



641og^(r™„) 5121og^(r: 



+ 



1536 



Notice that 



lim 

,in-*-0+ 



7{DP) 



C(3) 35 



32 768 



128 



0.006030910507 



+ . . <83) 



(84) 



is the corresponding Dirichlet sum rule for a unit circle. Once again we see that 
the derivative of Z^^^\rmin) diverges at rmin = 0. 
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The last case that we examine is sum rule of order four: for rmin 0+ it 
behaves as 



64 110592 1152 11520/ 

1 5 2147 677 



2561og^(r™„) 5121og3(r™„) 221184 log2(r„i„) 147456 log(r™„) 
2 { ^ 5 23 149 \ 

U4 V(r„„) ^ 2561og3(r„,„) ^ 3456 log2(r™„) ^ 138240 j ^ 



We have 



,Ji^o+^^''''^^'"™"^ ~ 1^ ^64 110592 +1152 + 115207 

w 0.0009438572210 , (86) 

which is the corresponding sum rule for a unit circle with Dirichlet be. The 
derivative of z[^^^ {rmin) diverges at r = 0. 

We now discuss the case of Neumann boundary conditions: the evaluation 
of these sum rules requires using the same expressions used for the case od 
Dirichlet be, with gH^^liXjX') instead of gi!^^l{x,x'). 

In the case of Neumann be and of a circular annulus with an infinitesimal 
hole we find 

2 [rminj ~ \^ 3g + g y V 48 96/ 



C(3) 3491 57r2 tt^ \ 



+ 



7 25 



321og2(r™„) 641og(r„i„)^ 

The logarithmic divergence of this expression signals the presence of eigen- 
values of infinitesimal magnitude as rmin — ^ 0+ . 

It interesting to study the sum rules of a circular annulus with an infinitesi- 
mal hole and mixed boundary conditions, either Neumann-Dirichlet or Dirichlet- 
Neumann on the internal and external borders respectively. 

For the Neumann-Dirichlet case we find the sum rule of order two: 

^2 \rn,rn) ~ \^ 43 33 j 43 + ^'"i" 48 288 

(3 11 
-log(r„i„) + — 

6377r^ 

- ^^-^L„log(r™in)(log(r„,„) + 4) + ... (88) 

We notice that Z^^^^ (rmin) tends to the corresponding sum rule for the 
unit circle with Dirichlet boundary conditions and that the first three derivatives 
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of Z^^^^ (jmin) arc finite at rmin = 0"*" and in particular that lim^ 
0: however the presence of a term r^„j„ log(r„ii„) implies that lim^^ 

00. 

Similarly we find the sum rule of order three: 



0+ dZf^^^{rrain)ldr, 

(r. 



>0+ " ^2 



nitn 
4 



Z. 



{NDP) 



iXmin ) 



C(3) 
32 



1781r"* 



128 ^ 768 
19 



7lH 



1536 



23040 



^''mmlog(r™i„) + r 



6 

rain 



7C(3) _ 197r2 
~32 



128 



( ^ log^(r™„) + ^ log(r,„i„) + ^ 



mm)l dr^^^ 



(89) 

162319\ 
92160 j 

(90) 



In this case we observe that Z^^^^^ (r„ 



^3 v'""") tends to the corresponding sum 
rule for the unit circle with Dirichlet boundary conditions and that the first 
three derivatives of Z^^^^\rmin) are finite at r,„i„ — 0+ and in particular 
that 1110.^^.^-^0+ dZ^^^^\rmin)/drmin — 0- however the presence of a term 
log{rmin) implies that limr„^„^o+ d'^Z^^^^\rm,in)/dr^in = oo. 



For the sum rule of order four we have 



Z, 



{NDP) 



03) 
64 



11520 1152 



3491 \ 
110592 j 



1691rL„ 
138240 



(91) 



40489r^_ 

829440 
13140797rL, 

116121600 



109rl„log( 



2304 



96 



,log^(?-« 



N , 571r^,„ log(r „^„) 
^' 1152 



(92) 



Once again we observe that Z^l^^^"* (r„ 



i) tends to the corresponding sum 
rule for the unit circle with Dirichlet boundary conditions and that the first 
three derivatives of z[^''^'^\rmin) are finite at Vmin — 0+ and in particular 
that liinr^^^^o+ dzl^''^'^\rmin)/drrnin — 0: however the presence of a term 
r^j„ log(r,„i„) implies that limr^^^^o+ d'^Z'-^'^^^i) /drf^^^ = oo. 

A different behavior is observed for the sum rules of the Dirichlet-Neumann 
case: for example, for the sum rule of order two we have 



y(DNP)( 



31og(r™i„) 



57r2 

Is" 



19 
32 



85rd 



48 



^(93) 



which diverges for 
boundary conditions. 



0+ , as in the analogous case of Neumann- Neumann 



5.2. Circular sector 

We consider the circular sector represented in FigHJ this domain is obtained 
applying the map ([76]) to the rectangle log r, — i log r] x [— 0], where r — >■ 
0+. Therefore it is natural to extend the analysis that we have done for the case 
of the circular annulus to this case. 
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The eigenfunctions of the negative laplacian on this domain are 

riTT , 



^n,ki'^,0) = N„k J !^{ankr) sin 



where Nnk is a normalization constant: 



V2 



(94) 



(95) 



The eigenvalues of the negative domain are then: 

Enk = «L , (96) 

where a„fc is the fc*^ zero of Ji^{x). 

We first discuss the case of t)irichlet boundary conditions at the borders of 
the circular sector; the general rules that we have derived in Section [3] can be 
applied straightforwardly. 

The sum rule of order two is 



oo 



n=l 
^2 



where V'*-"-' (z) = (— l)"+"'^n! X^fcLo ffc+zV+t ^ polygamma function. 

IS 

1 



^k=0 (fe+z)i 

The sum rule of order three is 



(7 + 2)' (^+4) (^+6) 



^ ^(^) + 1 



7<?!' (0) 2</> 



1287r 



The sum rule of order four is 



1287r \ TT 



(98) 
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967r^ ^ \ TT 
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5127r ^ V TT 
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127(/) 



20 



3847r2 
150 



46087r V TT / 5127r 



46087r 



(0) 



20^ 





TT 




-) 






(99) 
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Figure 2: Circular sector obtained acting with the map I I76II the rectangle logr, — i logr] X 
[-<!>, 4>], with r -i> 0+. 



Table 1: Dirichlct sum rules for circular sector at specific angular values 







TT 

4 

TT 

2 

371- 

4 

TT 


1 tt" log(4) 
32 ' 128 32 

TT^ 3 

96 ^ 32 
35 1 977" 31og(4) 
64 ^ 128 32 

TT^ 37 






TT 

4 

TT 

2 

3tt 
4 

TT 


7C(3) 7 37r" 1 log(4) 
512 768 1024 ' 64 

C(3) , 31 TT^ 

, , 64 1536 256 
189C(3) 6653 277r^ , 31og(4) 
512 26880 1024 J 64 
C(3) , 43 7r=^ 
8 ^ 7680 64 






TT 

4 
7r 
2 

371- 
4 

TT 


7C(3) 1 1 1 37r^ 1 77* 171og(4) 
1024 ' 36864 ' 2048 ' 24576 2304 
C(3) 1795 1 5t? I TT* 
128 110592 ' 2304 ' 23040 
189C(3) 256171 , 277r^ , 277r^ 171og(4) 
1024 1290240 ' 2048 ' 8192 768 
C(3) 33569 I 571-^ , 7r* 
16 430080 576 1440 
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5.3. Circular annulus with inhomogenous density 

We now discuss the case of a circular annulus with a density which depends 
only on the radial coordinate. For simplicity we assume an annulus of internal 
radius r„i„ and external radius r„iax = 1 a-nd with density 



ib + iU.rl^„-l)r'> 
2{r 



P(r) - ' \:\-' . (100) 



These drums all have the same mass M = 2n p{r)rdr — 7r(l — r^j„). 

The eigenmodes of these annular membranes are the eigensolutions of the 
Helmholtz equation 

(-A) e) = E,,p{r)i,r.{r, 9) (101) 

where r £ {r^in, 1) and 9 e (0, 27r). 

Using the conformal map ([7S)) we convert the Helmholtz equation (|10ip into 



d d 



$„(x,2/) = £;„E(x)$„(a;,y) (102) 



where x S (log(r„„;„)/2, - log(r™i„)/2) and y € (-7r,7r) and 

E(x) = (r™„e2-)p(V^e") (103) 
It is easy to check that 



E(x) = I](-x) 

b=-2+5 



b=-2-(5 



and therefore we conclude that the annuli with exponents b ~ —2±5 are isospec- 
tral. Using 6 = 2 we see that the uniform annulus is isospectral to the annulus 
with density p{r) = r'^^^Jr'^ (this case has been studied by Gottlieb in ref. (20jf^. 
The sum rule of order two for this inhomogeneous annulus is 



8(& + 2)M-™+'-l) V(-m») 
• [8 (r'+^n -lf + {b + 2) log(r™„) {-5r'^+* + (fe + 2) (^+4 + l) log(r™„) + 5) 

where we have defined 

^ {rLn if (-32n^ {n' - (b + 2f) ^ + {{b + 2f An'f r^ + {{b + 2f - 4^^)^) 



A more general class of radially isospectral annular membrane is discussed in ref. |2ll | 



23 



0.04 r 




-10 -5 5 

Figure 3: Sum rule of order 2 for a circular annulus with density p(r) 
The curves from top to bottom correspond to v min — 

0.1,0.25,0.5. 



b 

{b+2)(r 



2(r^+/„-l) 



+ \{h + 2) (d„. - 1)' r-?: {{h + n + 2){b + 2n + 2)^,^+^ + _ 2n + 2)^(6 - n + 2)) 
+ \ib + 2) {rl^^ - ifrll,, {{b -2n + 2f{b-n + 2)r,t+4 + (6 + n + 2)(6 + 2n + 2)^) (105) 
P„ ^ 2{b^ + Ab-An^ + A)\b' + Ab-n^ + A){ii+l-lf{r-^l^^r^:^,^:)' (106) 
In particular for b = —2 and rm^n — > 0^ wc find: 

360 81og(r™„) 3601og'(r„„;„) ^ ' 

In Figl3]we plot the sum rule Z2^\b) at three different values of rmin as a 
functions of 6. The curve is symmetric with respect to the axis 6 = —2, as a 
result of the isospectrality discussed earlier. 

In FiglHwe plot the sum rule z'^\—2) as a function of Tmin (solid line); 
the dashed line is the asymptotic behavior of eq. (|107p . 



6. Conclusions 



We have extended the results of Ref. [15[ to two dimensional domains of 
arbitrary shape and density; we have derived a general expression for the sum 
rule of integer order n in terms of the Green's functions of the homogeneous 
problem and we have proved that the sum rule of any integer order can be 
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Figure 4: (—2) as a function of r^i„ for a circular annulus with density p(r) = 

^ J J^^'" — The solid line is the exact result, whereas the dashed line is the asymptotic 
formula JTOTj. 

expressed in terms of a single series for the case in which the inhoniogeneity is 
only along a direction. We have also discussed the generalization of these results 
to higher dimensions. 

As an application of the formulas that we have derived, we have calculated 
the sum rules for an homogeneous circular annulus, for a homogeneous circular 
sector and for a circular annulus with a radially varying density. We have found 
explicit expressions for the sum rules of a circular annulus subject to different 
boundary conditions and with an infinitesimal internal hole: in the case of 
Dirichlet boundary conditions on the external border we have proved that the 
sum rules tend to the corresponding sum rules for the unit circle. 

The possibility of obtaining exact expressions for the sum rules of the eigen- 
values of the negative laplacian on arbitrary domains may be exploited to esti- 
mate the higher order coefficients in Weyl's asymptotic law. This aspect would 
provide a natural extension of the present work and we plan to consider it in 
the future. As we have discussed in [l^, following [sf, the sum rules may be 
used to provide rigorous upper and lower bounds to the lowest eigenvalue. 
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